Transport measurements across a tunable potential barrier in graphene 
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The peculiar nature of electron scattering in graphene is among many exciting theoretical pre- 
dictions for the physical properties of this material. To investigate electron scattering properties in 
a graphene plane, we have created a gate-tunable potential barrier within a single-layer graphene 
sheet. We report measurements of electrical transport across this structure as the tunable barrier 
potential is swept through a range of heights. When the barrier is sufficiently strong to form a 
bipolar junctions (npn or pnp) within the graphene sheet, the resistance across the barrier sharply 
increases. We compare these results to predictions for both diffusive and ballistic transport, as the 
barrier rises on a length scale comparable to the mean free path. Finally, we show how a magnetic 
field modifies transport across the barrier. 



The recent discovery by Novoselov et al. [ij of a new 
two-dimensional carbon material - graphene - has trig- 
gered an intense research effort 0, 0, Carriers in 
graphene have two unusual characteristics: they exist on 
two equivalent but separate sublattices, and they have 
a nearly linear dispersion relation. Together these in- 
gredients should give rise to remarkable transport prop- 
erties including an unusual suppression of backscatter- 
ing Unlike in conventional metals and semiconduc- 
tors, electrons in graphene normally incident on a po- 
tential barrier should be perfectly transmitted, by anal- 
ogy to the Klein paradox of relativistic quantum mechan- 
ics 0, 0, S] ■ Backscattering would require either breaking 
of the "pseudospin" symmetry between electrons living 
on the two atomic sublattices of the graphene sheet Q or 
a momentum transfer of order the Fermi wavevector kp, 
which can only be produced by a sharp, atomic-scale step 
in the potential. In addition to the obvious relevance to 
future graphene-based electronics, understanding trans- 
port across potential barriers in graphene is essential for 
explaining transport in graphene close to zero average 
density, where local potential fluctuations produce pud- 
dles of n- and p-type carriers Q. We report an experi- 
ment in which a tunable potential barrier has been fabri- 
cated in graphene, and we present measurements of the 
resistance across the barrier as a function of the barrier 
height and the bulk carrier density. 

To create a tunable potential barrier in graphene, we 
implemented a design with two electrostatic gates, a 
global back gate and a local top gate (Fig.[T^.) A voltage 
Vb applied to the back gate tunes the carrier density in 
the bulk of the graphene sheet, whereas a voltage Vt ap- 
plied to the top gate tunes the density only in the narrow 
strip below the gate. These gates define two areas in the 
graphene sheet whose densities n2 - underneath the top 
gate - and ni - everywhere else - can be controlled inde- 
pendently (Fig. [lb.) The graphene was deposited by suc- 
cessive mechanical exfoliation of natural graphite crystals 
using an adhesive tape (Nitto Denko Corp.) After ex- 
foliation, thin graphite flakes were transferred onto a chip 
with 280 nm thermal oxide on top of an n++ Si substrate, 
used as the back gate. A single-layer graphene sheet was 
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FIG. 1: a) cross-section view of the top gate device, b) sim- 
plified model for the electrochemical potential U of electrons 
in graphene along the cross-section of a). The potential is 
shifted in region 2 by the top gate voltage and shifted in 
both regions 1 and 2 by the back gate voltage, c) Optical 
image of the device. The barely visible graphene is outlined 
with a dashed line and the PMMA layer appears as a blue 
shadow. A schematic of the four-terminal measurement setup 
used throughout the paper is shown. 



identified using an optical microscope, and 30 nm thick 
Ti/Au leads were evaporated onto the sheet using stan- 
dard e-beam lithography. To form the dielectric layer 
for the top gate, a thin layer of PMMA (molecular mass 
950K, 2% in anisole) was then spun onto the Si chip at 
4000 rpm for 35 s and baked at 180°C for 2 minutes. 
On one section of the graphene sheet, the PMMA was 
cross-linked by exposure to 30 keV electrons at a dose 
of 21, 000 AiC cm"2 The unexposed PMMA was re- 
moved by soaking the chip for 10 minutes in acetone. 
Finally, a 50 nm thick, 300 nm wide Ti/Au strip was de- 
posited on top of the cross-linked PMMA using standard 
e-beam lithography to form the top gate (Fig. [TJ:.) Ra- 
man spectroscopy of the sheet indicates that it is in fact 
composed of a sing le layer [ll|,[il[ll. 

We performed conventional four-terminal current- 
biased lockin measurements (100 nA at 12.5 Hz) of the re- 
sistance of a PMMA-covered L ~ 1.3 /im by w = 1.7 /im 
section of a graphene fiake, in which a potential barrier 
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FIG. 2: a) Resistance across the graphene sample at 4 K as 
a function of the top gate voltage for several back gate volt- 
ages, each denoted by a different color, b) Two-dimensional 
greyscale plot of the same resistance as a function of both 
gate voltages. Traces in a) are cuts along the correspondingly- 
colored lines. 



can be induced by the top gate (Fig. [TJ;). All the mea- 
surements were performed in liquid Helium at 4 K. In 
Fig. [^b, this resistance is plotted as a function of both 
the back gate voltage H and the top gate voltage Vt [11] . 
Two clear white lines appear, indicating local maxima 
in the resistance in the different regions of the graphene 
sheet. In the bulk of the sheet (region 1), the average 
carrier density is given by ni = Cb{Vb — V^)/e, where 
the back gate capacitance per area Cb ~ 14 nF cm~^ is 
inferred from Hall measurements of other graphene flakes 
on the same oxide layer. The horizontal white line marks 
the neutrality point (ni = 0) in this region, allowing us 
to estimate « 31.5 V. In region 2, the average density 
n2 — ni + Ct{Vt — V°)/e is modulated not only by the 
back gate, but also by the top gate, with capacitive cou- 
pling Ct- The voltage is not necessarily zero, since 
the chemical doping in regions 1 and 2 can differ, and 
this difference can vary widely from sample to sample. 
The diagonal white line marks the neutrality point un- 
derneath the top gate (n2 = 0). The slope of this line 
provides the relative coupling of the graphene sheet to the 
two gates: Ct/Cg « 6.8, so Ct « 1.0 x 10^ nF.cm"^ Us- 
ing the PMMA thickness of 40 nm as measured by AFM, 
this value leads to a dielectric constant epMMA = 4.5 
(close to the accepted room temperature value for non- 
cross- linked PMMA). The crossing point of these white 



lines yields = -1.4 V [13| 



Transport fluctuations seen in Fig.[2]are a reproducible 
function of gate voltages and magnetic field (univer- 
sal conductance fluctuations [l3|) with amplitude SG ~ 
0.2e^/h [1^. The magnetoresistance is almost perfectly 
symmetric in magnetic field, as expected [13] • We ex- 
tract the phase coherence length « 4 /zm and the 
inter-valley scattering length Li_v ~ 0.15 /Ltm from the 
weak localization peak. Li^ > Li_v indicates that the 
sample is lying flat on the substrate [13, [H, [3. 

Fig. [^h- shows selected cuts through Fig. [^h: device re- 



sistance across the potential barrier as a function of the 
top gate voltage Vt, for several values of the bulk den- 
sity ni. Each curve has a maximum close to 712 — 0, 
arising from the enhanced resistivity of the graphene in 
region 2. However, each curve is noticeably asymmet- 
ric with respect to this maximum: the resistance de- 
pends on whether or not the carriers in region 2 (elec- 
trons or holes) are the same type as those in region 1. 
Specifically, for given absolute values of the densities 
Til and n2, the resistance is always higher if the car- 
rier types in the two regions are opposite (nin2 < 0) 
than if the carriers are the same throughout. In or- 
der to highlight the effects of pn junctions between re- 
gions 1 and 2, we extract for each value of ni the part 
of the resistance which depends on the sign of the 

carriers in region 2: i?„,(n2) = i?^7"(n2) + (^2), 
^nr(-^A = KTi^2) and (-712) = -i?°f (n2) 
(Fig. [3]) [20| . The resistance of the device away from the 
junctions (region 1 and possibly the interior of region 2) 
does not depend on the sign of n2 and hence is entirely 
contained in which we do not examine further. 
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FIG. 3: a) Odd component of the resistance: the part which 
depends on the sign of the density 712 in region 2. b) Greyscale 
plot of the odd component of the resistance for many values 
of ni . Colored lines are the cuts shown in a) . 



The presence of pn junctions between regions 1 and 2 
is associated with a substantial increase of the overall re- 
sistance (white areas compared to black areas of Fig.ISj)). 
We label each section of Fig. [5)3 with the carrier types 
in region 1 and 2, "p" for holes and "n" for electrons, 
to emphasize that enhanced resistance is associated with 
an npn or pnp junction. The junction-sensitive resis- 
tance curves are almost symmetric upon simultaneous 
sign change of both densities ni and 712: left-right re- 
flection and color swap in Fig. [3^, or 180 degree rotation 
about the center in Fig.[3)D. Deviations from this symme- 
try are presumably associated with uncontrolled spatial 
fluctuations in the density. Strikingly, a sharp step ap- 
pears in the resistance as the boundary between nnn and 
npn or pnp configurations is crossed. We can explore the 
underlying physics of such a resistance increase in two 
opposite regimes: strongly diffusive or ballistic. Which 
applies to experiments depends on whether the elastic 
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mean free path is greater than or less than the length 
over which the potential barrier rises. 

In a simple model, the electrostatic potential U{x) 
induced by the set of two gates in the graphene sheet 
changes linearly between regions 1 and 2 over a width 2d 
(Fig. [Dd). Assuming perfect screening in the graphene 
sheet (a good approximation for our relevant density 
range), we estimate d « 40 nm independent of carrier 
density, by solving the Laplace potential in the region 
between the top gate and the graphene. If the elastic 
mean free path 1^ is much shorter than d, the total re- 
sistance can be estimated by integrating the local resis- 
tivity across the device. The conductivity cr(x) at any 
position X depends only on the density of charge carriers 
n{x) and can be inferred from resistance measurements 
where the density is uniform [l3l |. It can be approxi- 
mated by an interpolation between low and high density 

behavior (7(2;) w {efin(x)) +0-^111 where the mo- 

bility is /i « 2 X lO'^cm^V^^s^^ and fJmin ~ 4e^/ft-. Using 
d = 40 nm and a width w = 1.7 /im, we plot the predicted 



odd part of R{Vt) with respect to the voltage 
each value of ni on Fig. [4^. 
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FIG. 4: a) Using the strongly diffusive model described in 
the text, one can predict the resistance as a function of the 
top gate voltage Vt for several values of the density ni (each 
represented by the same color as in Fig. |3]). Here we plot 
the odd part of this calculated resistance (cf. Fig. |3} for a 
barrier smoothness d = 40 nm and a width w = 1.7 ^m. b) 
Similar curves for several densities n\ using the ballistic model 
described in the text. The curves are plotted only for densities 
n2 not too close to zero: 712 > 10^^ cm^^ (the model diverges 
at 712 = 0) and dashed lines link the curves at opposite sides 
of 712 ~ 0. 



At high enough density, where the conductivity a of 
graphene is proportional to the density of charges, one 
can define a mean free path 1^ in an isotropic diffusive 
model by cr = 2fci?/ee^//i, where kp — nn is the effective 
Fermi wavevector defined relative to the K-point of the 



Brillouin zone [|21]. For m 



n2 
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see in Fig. that a « 6^- [221 . Therefore, le « 0.03 ^m 



for n 



1 X 10 cm . The mean free path is then 



gate. Hence, the diffusive model is a poor approximation 
and is unable to reproduce the data shown in Fig. [S] 

The opposite regime of ballistic transport has been 
considered in the limits of sharp and smooth potential 
steps. We use here the results of Cheianov et al. [1], 
which are valid in the limit |nij, \n2\ ^ d^^. Note that 
the results for a sharp barrier [6| lead to the same qual- 
itative results for R'^^'^{n2) but are more than 10 times 
smaller than what follows. The transmission probability 
t{9i) for electrons impinging from region 1 on the in- 
terface between regions 1 and 2 depends strongly on the 
angle of incidence 9i with respect to normal incidence, 
and is given by 



if 7iirt2 < 0. (1) 



We incorporate both interfaces of the potential barrier 
into the resistance calculation in the following way. In 
the fully ballistic regime {le ^ L2,d, where L2 is the dis- 
tance between the two interfaces), all possible reflections 
at the interfaces, and their interferences, must be taken 
into account. However, in the regime of this experiment, 
where le ^ d <^ L2, the two interfaces can be considered 
as independent resistors in series. Using a Landauer pic- 
ture, the (antisymmetric) conductance through the full 
barrier is then 



mi, — — m" 



..l|l/2 + |„2|l/2 ^ 



(2) 



already of the order of d for a few volts applied to the back 



where we sum the transmission coefficients of inci- 
dent modes having all possible values of the quan- 
tized wavevector component perpendicular to the in- 
terface ky — 2T:my/w, with my integer and m^^^ = 

■\/7rmin(|7T,i I, \n2\)w/2TT . Each such mode carries a con- 
ductance 4e^//i, where the degeneracy factor 4 is charac- 
teristic for graphene. The presence of two interfaces in 
series gives an overall prefactor of 1/2. 

Using d = 40 nm and w — 1.7 /xm, we plot the calcu- 
lated odd part of the resistance in our ballistic model as 
a function of Vt in Fig. IHd for the same densities ni as in 
the strongly diffusive case, assuming that the resistance 
in the case where the sign of n is uniform is much smaller 
than l/Gnpn- Although the results of the ballistic model 
are of the same order of magnitude as the measured val- 
ues of the odd part of the resistance, the model is unable 
to reproduce some aspects of the experimental data. In 
particular, the model does not explain why R^^'^ does 
not only jump but continues to increase as \n2\ passes 
beyond zero density. This behavior is particularly sur- 
prising as, for fixed n2, R°^'^ decreases with increasing 
|7ii|. Furthermore, by construction, the model does not 
apply close to zero density ni or ?i2. 

A complementary test of the unusual transmission 
through a potential barrier in graphene as a function of 
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the angle of incidence was proposed by Cheianov et al. . 
In the fully ballistic regime, the resistance across the bar- 
rier is predicted to increase as soon as a magnetic field 
applied perpendicular to the graphene sheet gets higher 

/ \l/2 

than a typical value w h/e I y^TT\n2\/TTL2d\ , where 

L2 is the barrier length (0.3 /im via SEM). Using the val- 
ues d — iO nm, L2 = 0.3 jum, and 712 = 2 x lO^^cm^^, 
one gets « 0.1 T. 
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height is not yet understood. A clear explanation of this 
behavior is essential for realizing many exciting proposed 
graphene electronics applications, such as electron focus- 
ing [13]. 
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FIG. 5: Odd part of the resistance as a function of the back 
gate voltage Vb, with the constraint n2 = — ni. Measurements 
were taken at 4 K for several values of the perpendicular mag- 
netic field B. 

Fig. [5] shows measurements of the antisymmetrized re- 
sistance across the device for several values of a magnetic 
field B applied perpendicular to the graphene sheet, with 
the top gate voltage Vt adjusted to maintain n2 ~ —ni. 
This odd part of the resistance shows an enhancement 
as the magnetic field increases above 2.5 T, an order 
of magnitude higher than the predicted transition field 
B^,. The fully ballistic model is therefore unable to ex- 
plain the results of our measurements in magnetic field. 
In fact, the diffusive model also predicts an increase of 
Rn'l'^ as B increases, due to magnetoresistance of the 
strip near zero density between p and n regions. If 
the entire graphene sheet is set to zero density, its re- 
sistance increases with magnetic field similarly to -R°'J'^ 
at 712 ~ -ni w 2 X lO^i^cm^^ (pjg ^ ^phis qual- 
itative agreement indicates that while ballistic physics 
plays a role in transport across the individual pn junc- 
tions, transport through the full npn barrier is far from 
ballistic. 

In conclusion, we have fabricated a gate-tunable bar- 
rier device from a single-layer graphene sheet, and have 
hence created bipolar junctions within graphene. We 
study transport across the potential barrier as a function 
of Fermi level and barrier height, as well as in the pres- 
ence of an external magnetic field, and demonstrate that 
a sharp increase in the resistance occurs as the potential 
crosses the Fermi level. This increase is better described 
by a ballistic than a diffusive model, but the dependence 
of the additional resistance as a function of the barrier 
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